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Abstract 



>, 

A light sterile neutrino species has been introduced to explain simultaneously 

the solar and atmospheric neutrino puzzles and the results of the LSND ex- 

OS . periment, while providing for a hot component of dark matter. Employing 

this scheme of neutrino masses and mixings, we show how matter-enhanced 

q . active-sterile {v^ T v^ v s ) neutrino transformation followed by active-active 

\ v u.r ^ v e) neutrino transformation can solve robustly the neutron deficit 

problem encountered by models of r-process nucleosynthesis associated with 

neutrino-heated supernova ejecta. 
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I. INTRODUCTION 

In this paper we show how the invocation of sterile neutrinos in a novel transformation 
scenario can enable the production of the heavy rapid-neutron-capture (r-process |T|) ele- 
ments in neutrino-heated ejecta from supernovae. Interestingly, the current hints of (and 
constraints on) neutrino oscillation phenomena from considerations of solar || and atmo- 
spheric neutrinos || and from the LSND experiment 0] are difficult to explain with only 
two independent neutrino mass-squared differences (corresponding to three active neutrino 
species) [^|. By contrast, the data are explained readily in terms of neutrino oscillations 
with three independent mass-squared differences, corresponding to four neutrino masses 
f||7j. However, the observed width of the Z° || can accomodate only three light, active 
neutrinos. Therefore, a fit to the data requires introducing at least one light "sterile" neu- 
trino species which does not have normal weak interactions [e.g., a Majorana SU(2) singlet 
neutrino] . 

The most promising site for r-process nucleosynthesis is the neutrino-heated material 
ejected relatively long (~ 10 s) after the explosion of a Type II or Type Ib/c supernova 
||. However, detailed calculations of the conditions which obtain above the neutron star 
remnant in such r-process models show that the neutron-to-seed ratio, R, is too low to allow 
the production of the heaviest r-process species ITDI. (The "seed" nuclei which capture 



neutrons to make the heavier species have nuclear mass numbers between 50 and 100.) We 
require R > 100 to effect a good r-process yield for heavy nuclear species, but the models 
with conventional neutrino physics and conventional equations of state for nuclear matter 
all give smaller values of R. This is the neutron deficit problem for these models of the 
r-process. 

In all models, R is determined by the net electron-to-baryon number (i.e., the electron 

fraction), Y e , the entropy-per-baryon in the ejecta, S, and the dynamic expansion timescale, 

tdyn, associated with the material ejection process [|10j- Though general relativistic effects 

]"T| and multi-dimensional hydrodynamic outflow |12[] have been invoked to increase 5* and 



decrease tdyn (both of these changes favoring larger R) enough to solve the neutron deficit 
problem, these solutions are at best finely tuned. On general grounds we could argue that 
the only robust way to obtain R > 100 in conventional neutrino-heated ejecta is to decrease 
Y e and/or maintain its low value by invoking new neutrino physics, e.g., introducing a sterile 
neutrino as discussed in this paper. (By "robust" we mean robust to astrophysical uncer- 
tainties in the detailed characteristics of neutrino- heated outflow.) We cannot definitively 
proclaim at this point that new neutrino physics is required to understand the production of 
heavy r-process elements. Such a proclamation calls for the accomplishment of the following. 
First, we would have to establish that at least some of the r-process material originates in 
conventional neutrino-heated ejecta. Second, we would have to understand the thermal and 
hydrodynamic evolution of the very late stage neutrino-driven "winds" from a proto-neutron 
star. 

The proposition that some of the r-process material comes from environments with in- 
tense neutrino fluxes is supported by recent studies of neutrino effects during and immedi- 
ately following the r-process [|T3|JH|] . Perhaps more to the point, recent observational data on 
the abundances of r-process species in old, very metal-poor halo stars in the Galaxy seem 
to be quite consistent with what is expected from the r-process scenario associated with 



neutrino- heated ejecta ||15|| . However, the second issue we would have to resolve is a vexing 

one. Despite the extensive numerical insights into the early phase of the supernova evolution 

16| , consistent and accurate hydrodynamic scenarios for neutron stars to lose mass through 



neutrino heating at very late times remain to be established. On the positive side, Burrows 



1 17 1 has recently followed his supernova calculations long enough to see the formation of a 
neutrino-driven wind at late times. In the region relevant for nucleosynthesis, the wind in 
his calculations resembles in broad brush the outflow models studied here. 

This paper is organized as follows. In Sec. II, we describe a simple model for the 
neutrino-heated outflow and discuss the associated alpha effect which causes the neutron 
deficit problem for r-process nucleosynthesis. In Sec. Ill, we discuss the treatment of 
active-sterile plus active-active neutrino transformation in supernovae and describe how 



such transformation can evade the alpha effect. The consequence of such transformation for 
the evolution of the electron fraction Y e in the neutrino-heated outflow is studied in Sec. IV. 
Conclusions are given in Sec. V. 

II. NEUTRINO-HEATED OUTFLOW AND THE ALPHA EFFECT 

In what follows we adopt a simple exponential wind model for the neutrino-driven outflow 
above the surface of the hot proto-neutron star produced by a supernova explosion. In this 
model the enthalpy per baryon is roughly the gravitational binding energy of a baryon, 
leading to a relation between radius r& (in units of 10 6 cm), entropy-per-baryon S^oo (in 
units of 100 times the Boltzmann constant), and temperature Tg (in units of 10 9 K): 

^ J2^_ / M NS \ 
r6 ~T 9 S W0 {lAM Q )> [L) 

where Mns is the mass of the proto-neutron star. Therefore, for a constant entropy per 
baryon characterizing the adiabatic expansion of the outflow, the temperature parameterizes 
the radius |11|| . The entropy typically will be carried almost exclusively by relativistic 
particles (photons and electron-positron pairs) so that S ~ (2ii 2 /4:5)g s T 3 /(pNA), where g s 
is the statistical weight for relativistic particles and N& is Avogadro's number. This leads 
to a simple relation between (matter) density p and radius: 

p « 3.34 x lOVm- 3 ^)^^ « 3.8 x ^\c^(J^J (^S^rf, (2) 

where we have scaled g s assuming its value at T g > 10. As we will be mostly interested in 
processes occurring at T 9 > 10, the dependence on g s will be suppressed hereafter. 

In the exponential wind the radius of an outflowing mass element is related to time t by 
r = r exp [(£ — £o)/tdyn]> Here tdyn is an assumed constant material expansion timescale. 
This implies an outflow velocity proportional to radius, v = t/tdyn, so that the neutrino- 
driven wind will remain roughly self similar for timescales on which the neutrino luminosities 
and energy spectra and the neutron star radius can be regarded as constant. Clearly, at 



some point above the neutron star surface the exponential wind must go over to a linear 
expansion of the radius, i.e., a "coasting" outflow. Nevertheless, the exponential wind regime 
should encompass the region above the neutron star where most of the biggest obstacle to 



successful r-process nucleosynthesis, the alpha effect (see below), is operative |T3 . 

In neutrino-heated ejecta, neutrino interactions with matter supply the requisite energy 
for ejection of nucleosynthesis products. The total amount of heating through these interac- 
tions determines S and tdyn in the ejecta. Most of this heating occurs close to the neutron 
star and S and tdyn are set at Tg ~ 20. However, the dominant interactions for neutrino 
heating 

v e + n -> p + e~ (3) 

and 

^e + V ~ * n + e+ (4) 

have prolonged effects on the neutron-to-proton ratio (n/p = 1/Y e — 1) in the ejecta. First 
of all, in the region where T 9 > 10 and free nucleons are favored by nuclear statistical 
equilibrium (NSE), the competition between the processes in Eqs. (^) and (|J) leads to n/p ~ 
Ap eP /A„ e n ~ (L Pe (E Pe )) J (L Ue (E Ue )). Here X PeP and X UeH are the rates for the reactions in Eqs. 
(|) and (H), respectively, L Pe and L Ue are the v e and u e energy luminosities, respectively, and 
{E Pe ) and (E Ve ) are the average energies characterizing the corresponding neutrino energy 
spectra. Absent neutrino oscillations (flavor/type mixings), we expect (E u ) w (E p ) m 
(E Ut ) p» (Ep T ) > (Ep e ) > (E„ e ) and, hence, n/p > 1 (Y e < 0.5, corresponding to neutron- 
rich ejecta) |18| . We note that the numerical value of the n/p ratio depends on careful 



calculation of the rates X PeP and A„ e „ |19[ and accurate determination of L Pe , L Ue , (E Pe ), and 



(E Ue ) IpOfl . However, the alpha effect discussed below is insensitive to the value of the n/p 
ratio obtained at T 9 > 10. 

The equilibrium of the n/p ratio with the v e and v e fluxes described above is maintained 
until the ejecta passes the weak freeze-out radius (at Tg ~ 10), rwFO, beyond which X Uen , 



Xp eP < 1/tdyn and free nucleons are no longer favored by NSE fl8fl . However, even beyond 
r WFO) v e capture on neutrons can force down the n/p ratio and, hence, the neutron-to-seed 
ratio R. This is the so-called "alpha effect" pl| , |13|| : as material flows to regions where the 



composition favored by NSE shifts from free nucleons to a mixture of free nucleons and 
alpha particles, protons will be incorporated into alpha particles, leaving a disproportionate 
fraction of neutrons exposed to the intense flux of v e . The process in Eq. (||]) then converts 
some of these neutrons into protons (which are then immediately incorporated into alpha 
particles), thus progressively lowering the n/p ratio and, ultimately, R. 

The cumulative damage done to the n/p ratio by v e capture on neutrons beyond the weak 
freeze-out radius r W FO can be estimated crudely by integrating dY e /dt « X Ven {t){\ — 2Y e ) to 
obtain 

Y e {t) « 0.5 + [F e (0) - 0.5] exp {-A„ en (0)r DYN [l - exp(-2t/r DYN )]}, (5) 

where initial values for the v e capture rate on neutrons and the electron fraction, X Uen (0) and 
Ye(0), respectively, can be approximated as their values at 7"\yFO Equilibrium of the n/p ratio 
with the u e and u e neutrino fluxes at early times requires that A^ en (0)r DY N ~ 1- Furthermore, 
the product A^ en (0)rDYN should be roughly constant for all models with conventional neutrino 
physics, since increased (decreased) neutrino luminosity will tend to increase (decrease) \ Vf , n 
and decrease (increase) tdyn proportionately. [Neutrino heating through the processes in 
Eqs. (HD and (|), which set Y e , is the principal determinant of t DY n-] Therefore, from Eq. 
(§) it can be seen that regardless of the freeze-out value Y e (0), sufficient exposure to the v e 
flux above r WF o will drive Y e close to 0.5 {n/p = 1, i.e., no neutron excess). 

In summary, a great paradox exists for r-process nucleosynthesis in a neutrino-driven 
wind. In such a scenario the neutrinos must supply enough energy [largely through the 
processes in Eqs. (^) and (|4])] to lift baryons out of the gravitational potential well of the 
neutron star. The gravitational binding energy per nucleon near the neutron star surface 
is ~ 100 MeV. Since the average v e and v e energies are ~ 10 MeV, a nucleon must interact 
with neutrinos ~ 10 times in order to acquire enough energy for ejection from the neutron 
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star. In turn, this requires a sufficiently large value of Aj, en (0)rDYN, thus providing conditions 
for a pernicious alpha effect which will cause a neutron deficit. Such a neutron deficit will 
preclude a successful r-process, especially the synthesis of the heaviest nuclear species. 

However, it is obvious that removal of the v e flux could neatly solve the neutron deficit 
problem by unbalancing the competition between the processes in Eqs. (^) and (^) in favor 
of neutron production (thus lowering Y e ) and by disabling the alpha effect in regions beyond 
twfo (thus maintaining a low Y e ). Furthermore, as neutrino heating is essentially completed 
at Tg ~ 20 (quite close to the neutron star surface), the beneficial effects of removing the v e 
flux can be obtained without affecting the general thermal and hydrodynamic characteristics 
of the neutrino-heated outflow (see Sec. III). 

III. EVASION OF THE ALPHA EFFECT VIA ACTIVE-STERILE PLUS 
ACTIVE- ACTIVE NEUTRINO TRANSFORMATION 

The existence of at least one light sterile neutrino could provide a means to reduce the 
v e flux at sufficiently large radius to leave the process of neutrino-heated ejection unhin- 
dered, yet sufficiently near or inside r WFO so as to disable the alpha effect, and thereby fix 
the neutron deficit problem for r-process nucleosynthesis. This r-process solution can be 
obtained by having an active-sterile (z/^ T ^ v s ) neutrino mass-level crossing followed by an 
active-active (z/^ T # v e ) neutrino mass-level crossing. (Hereafter, v s and v s denote left- and 
right-handed Majorana sterile neutrinos, respectively.) The first mass-level crossing converts 
z/ M and z/ T , which are emitted with the highest average energy from the neutron star, into 
harmless sterile neutrinos. Without this, significant conversion of v^ and v T into v e would 



drive the neutrino-heated supernova ejecta proton rich |L8[. With the supernova v^ and v T 
"sterilized" by the first mass-level crossing, the second mass-level crossing now only acts to 
convert v e emitted from the neutron star into v^ and v T . Charged-current capture reactions 
on neutrons are energetically forbidden for supernova v^ and v T . The net result is that the 
v e flux is reduced or removed, and the alpha effect cannot operate. Therefore, the neutron 



excess is preserved or possibly enhanced in this scenario. 

The sequence of neutrino mass-level crossings described above can occur in a four- 
neutrino scheme || which has two nearly degenerate neutrino doublets separated by a mass- 
squared difference chosen to be compatible with the LSND v^ v^ u e signal. In this scheme 
the lower mass neutrino pair would either give an active-sterile u e v^ v s mass-level crossing 
in the sun, or provide a "just so" vacuum mixing solution with 5m? es ~ 10~ 10 eV 2 to the solar 
neutrino problem ||22|| . For solving the solar neutrino problem via matter-enhanced v e ^ v s 



2 ^ in-5„Ar2 



mixing, the mass-squared splitting for the lower mass neutrino pair must be 5m es < 10 eV 



22fl . The mass-squared splitting for the higher mass neutrino pair is chosen to explain the 
observed atmospheric v^/v e data via maximal v^ v^ v T vacuum mixing. 

Figure 1 shows our adopted neutrino mass scheme. The mass-squared splittings shown 
are chosen to meet all experimental constraints and to solve the neutron deficit problem 
for r-process nucleosynthesis associated with neutrino-heated supernova ejecta. We would 
require m Vs > m Ue for matter-enhanced v e ^ v s mixing in the sun, but a "just so" solution 
to the solar neutrino problem could have m Vs < m Ve as well. (Issues of compatibility with 
Big Bang Nucleosynthesis limits aside, note that some kinds of singlet "sterile neutrinos" 
could possibly evade bounds on the "just so" u e ^ v s vacuum mixing solution stemming 



from the SuperK and the Chlorine Experiments p3|.) In our overall mass scheme z/ M and 



v T 



could share the role of providing a hot dark matter component [24 . 

Of the possible four-neutrino mass patterns, this one is most successful in evading con- 
straints. It has been shown that limits from accelerator and reactor experiments disfavor 



having one dominant neutrino mass, e.g., a3 + loral + 3 arrangement ||25|| . A two-doublet 
scheme with z/ M # v s mixing (the higher mass pair) explaining the atmospheric neutrino 
anomaly and v e ^ v T mixing (the lower mass pair) explaining the solar neutrino puzzle is 
in trouble with Big Bang Nucleosynthesis bounds |[26|j . 

With an appropriate choice of mixing angles and mass-squared difference for the splitting 
of the doublets, the neutrino mass scheme in Fig. 1 will lead to an efficient matter-enhanced 
^V,r ^ v s transition above the neutron star surface (and the neutrino sphere), yet below 



the region where an ordinary Mikheyev-Smirnov-Wolfenstein (MSW) |27]] matter-enhanced 



>,T 



v e mass-level crossing would occur [2c . Furthermore, the mass- level crossing for v e 



with energies most relevant for determining Y e can still lie in the region near or below the 
weak freeze-out radius rwFO- 



A. General Description of Neutrino Flavor/Type Evolution in Supernovae 

As discussed in Sec. IIIB, the overall problem of active-sterile plus active-active neutrino 
transformation in our case can be treated as that of two separate mass-level crossings, with 
each mass-level crossing involving effectively only two neutrino flavors/types. The neutrino 
flavor/type evolution through each mass-level crossing is governed by a Schroedinger-like 
equation: 



■A 

'dt 



a x {t)~ 


= (H v + H e + H vv ) 


'a x {t)~ 


a y (t)_ 




. a y(t). 



(6) 



where a x (t) and a y (t) are the time-dependent amplitudes for the neutrino to be in, for ex- 
ample, the flavor eigenstates \u e ) and |z^), respectively. The propagation Hamiltonian in 
Eq. (||) is given as the sum of three terms H v , H e , and H vv resulting from vacuum masses, 
forward scattering on electrons and nucleons, and forward scattering on "background" neu- 
trinos, respectively. 

In general, we can express the neutrino propagation Hamiltonian as: 

Y /-Acos29 + A + B Asin20 + 5 off 
H v + H e + H uv = - 

2 \ Asm26 + B* s Acos26-A-B 

where A = 5m 2 /(2E l/ ), with E v the neutrino energy and 5m? the appropriate vacuum 
neutrino mass-squared difference, A is the appropriate electron/nucleon background contri- 
bution, and B and B Q s are the diagonal and off-diagonal contributions, respectively, from 
the neutrino background. (The terms B g and B* s are complex conjugates of each other and 
vanish in the case of active-sterile neutrino transformation.) Note that a mass-level crossing 
or resonance will occur if the following condition is satisfied: 



(7) 



9 



Acos26 = A + B. (8) 

Table I gives appropriate expressions for the weak potentials A and B for all possible cases 
of 2 x 2 active-sterile and active-active neutrino mixings in terms of density p (in gem" 3 ), 
Avogadro's number Na, Y e , and the effective neutrino numbers for each of the three active 
species relative to baryons, Y Ve , Y v , Y Vt . 

The effective neutrino number for species a [a = e, //, r) relative to baryons is Y Va = 
(n® — n£ ff )/ pNa, where for example, the effective number density of neutrino species a at 



position r encountered by a "test" neutrino traveling in a pencil of directions is [18,|29,30 



1 L r°° r HO' - - 

< - E ^^ I <IE„, J - U (E,)P^(E„, r, tf)(l - n- • «). (9) 

Here /3 runs over all neutrino species to be considered (including species a), f Vp {E V0 ) is 
the normalized energy distribution function for neutrino species j3 at its birth position (the 
neutrino sphere of radius R u ), and P U/3 ^„ a (E Ul3 ,r,0') is the probability for an initial up 
to appear as a u a when encountering the test neutrino at position r. In general, this 
transformation probability will depend on the energy E Ug and direction O' of the background 
neutrino up, as well as on the position r. The last factor in Eq. (H) contains additional 
dependence on the direction O' of the background neutrino. 

For efficient conversion of neutrino flavors/types, evolution of the neutrino amplitudes 
through the mass-level crossing (resonance) must be adiabatic. In turn, this requires that 
the width of the resonance region be large compared with the local neutrino oscillation 
length. The width of the resonance region is 

5r = Ht&n29, (10) 

where Ti is the scale height of the weak potential at resonance: 

d\n(A + B) 



n 



(ii) 



dr 

The adiabaticity parameter 7 characterizing the evolution of the neutrino amplitudes 
through resonance is proportional to the ratio of the resonance width to the local oscil- 
lation length: 
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7 



Sm 2 sin 2 28 



H^4.6 



Sm 2 \ /25MeV\ /sin 2 2# N 



2E v cos28" ~\QeV 2 \ E u )\ 10" 3 ; V7.. r ) km 



n 



(12) 



Note that 7 > 3 corresponds to better than 99% conversion of neutrino flavors/types. 



B. The Case of Active-Sterile plus Active-Active Neutrino Transformation 



In the context of the general scheme in Fig. 1 for neutrino masses and mixings we 
make several specifications to facilitate our goal of removing the bulk of the v e flux in the 
appropriate region above the neutrino sphere. First we invoke maximal mixing between 
Up and v T , which is consistent with the SuperK atmospheric neutrino data. Furthermore, 
we assume that the vacuum mixing between v^ T and u e>a is small. In particular, we con- 
sider the following mixing scheme between the flavor/type eigenstates and the vacuum mass 
eigenstates: 

\ fWi)\ 
W2) 



( \u e ) \ ( cos 



\Wt)J 



(13) 



sm <p cos to sm <p sm to 
sin0 cos 4> cos u cos sin oj 
-sincj/v/2 coscj/v/2 l/V? \u 3 ) 

V sinuj/\/2 -COSLU/V2 I/V2J \\ui)J 

In Eq. (|13|), the vacuum mixing between u e and v a is largely governed by the angle 0, while 
that between v^ and v T is chosen to be maximal by our assumption. The angle uo essentially 
specifies the vacuum mixing between z/ MjT and z/ es . For small vacuum mixing between these 
two doublets, we require u -C 1. In Eq. (|T3"|), all the CP- violating phases are ignored. 



With the definition of 



1 *\ 



and 



\K) 



V2 ' 



V2 



(14) 



(15) 



Eq. fll3|) can be rewritten as 
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(W)\ 
if) 

11/ 

WO J 



( cos0 sin0cosc(j sin0sinc<j 0\ 

- sin cos cos u cos sin ui 

— sin uj cos cj 

V iy 



/h>\ 
V21 



(16) 



It is clear from Eq. ([TB]) that \u*) is a mass eigenstate in vacuum, i.e., H v \u*) oc |z/*). 
In fact, 1 1/*) is also an effective mass eigenstate in the presence of electron/nucleon and 
neutrino backgrounds, i.e., (H v + H e + H vv )\v*) oc \v*). This is because the amplitude of 
forward scattering on electrons and nucleons is the same for v^ and v T while the effective 
number densities of v^ and v T for neutrino-neutrino scattering above the neutrino sphere 
are the same due to the symmetry in transformations concerning v^ and v T . Therefore, the 
evolution of \u*) is decoupled from that of \v s ), \u e ), and |z/*). 

As discussed above, maximal mixing between v^ and v T allows us to reduce the problem 
of 4x4 neutrino mixing into that of mixing among u s , v e , and u* In the region above, yet not 
too far away from the neutrino sphere, this 3x3 neutrino mixing problem is characterized 
by the mass-level crossings for v* v^ v s and v* v^ v e transformations. Note that although 



8ml* s ~ 



8m 2 , in our neutrino mass scheme shown in Fig. 1, the weak potentials A and 



/i*e 



B in the neutrino propagation Hamiltonian are quite different for v* v^ v s and v* 



v„ 



transformations. The effect of the weak potential B from the neutrino background will be 
addressed in Sec. IV. For simplicity, we will neglect B in the following discussion. According 
to Eq. (H) and Table I, the mass-level crossing (resonance) for v* # v s transformation occurs 
when 



8m 



fl*S 



2E„ 



cos2# M * s = V^GppN; 



l-Y P 



while that for v* v^ v e transformation occurs when 



8m 2 ,, 



2E, 



¥-coa26u* e = V2G F pN A Y e . 



(17) 



(18) 



In Eqs. (|TTD and ([3~8|) , Gp is the Fermi constant, and O^g and # M * e are the appropriate 



two-neutrino vacuum mixing angles. From Eq. 



6,j,*e ~ usincj) <C 1. 



we have 6^* s ~ 



o;cos0 <C 1 and 
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In our proposed scheme to enable the r-process in neutrino-heated supernova ejecta by 
disabling the alpha effect, we require that the v* ^ v s conversion of relatively high energy 
neutrinos take place well below the weak freeze-out radius. The temperature at which weak 
freeze-out in the ejecta occurs is T 9 WFO « 10. According to Eqs. (|2|) and (|ITD , conversion 
below (at smaller radii than) the weak freeze-out radius then requires 



t wfo\3 Ev a(i_y e )/ 2 



K^BeV'^U^j 



SToo, (19) 



0.3 

where we have taken g s = 11/2, consistent with the conditions in the ejecta for Tg > T g WFO , 
and scaled the result assuming a typical value of Y e = 0.4 as obtained in numerical supernova 
models in the absence of neutrino transformation. Similarly, for the ^* v^ v e resonance to 
occur below the weak freeze-out radius, we require 

'T 9 WFO \ 3 / E v \(Y e 



,5m 2 >5.1eV 2 ^ _^_ JlfL £-£. (20) 

Me ~ V 10 / V25MeV/V0.4/ 10 ° v ; 

Therefore, a value of 8m 2 , s ~ 8m 2 , e ~ 6 eV 2 consistent with the LSND data would fulfill 
the requirements in Eqs. flTPD and (f2"U|). 

We note that for Y e > 1/3, the v* # v s resonance will occur before (at higher density 
and temperature than) the v* v^ v e resonance for a given neutrino energy [cf. Eqs. (|T^ ) and 
(0)]. Specifically, the temperature at resonance is 



8m 2 s \ l/ - i /25MeV\ 1 ^ 



»-~ u - 6 U# V-nr) 



0.3 



(i-n)/2 



1/3 



,1/3 



« (21) 



for the y* ^ i> s case, and 

— &) I/3 (^) 1/3 (^)^ 

for the i/* v^ v e case. According to Eq. ([I]), the radial separation between the two resonances 
has a typical value of 2 km. As shown below, this is much larger than the widths of both 
resonance regions. 

The scale height of the weak potential for both v* v=^ v s and v* v^ v e resonances can be 
approximated as 

13 



n 



<ilnp 



dr 



-l 



r -y, (23) 



where r res is the radius for the relevant resonance, and we have neglected the change in Y e 
compared with that in p. To ensure adiabatic conversion and accommodate the LSND data 
at the same time, we require sin 2 29 ~ 10~ 3 for 9 = 9 fl * s , 9^* e . [This is easily achieved by 
having uj ~ 10~ 2 and ~ 7r/4 in Eq. flTB|).] From Eq. (|TUD, the widths of the v* ^ v s and 
u* t=^ v e resonance regions are (5r)^* s ~ (<5r) M » e ~ 0.2 km, much smaller than the separation 
between the two resonances. Therefore, we can treat the evolution of v e as unaffected by 
the v* v^ v s resonance and that of v s unaffected by the u* v^ v e resonance. 

So far we have restricted our discussion to the case of Y e > 1/3. As the electron fraction 
Y e gets close to 1/3, the v* v^ v s and v* ^ v e resonances will approach each other and 
overlap. In addition, the weak potential governing the evolution of neutrino amplitudes will 
be dominated by the neutrino background for Y e w 1/3. However, as discussed in Sec. IV, 
the transformation of neutrino flavors/types, coupled with the expected rapid expansion of 
the neutrino-heated ejecta, will cause the actual value of Y e at a radius to differ significantly 
from the equilibrium value corresponding to the local v e and v e fluxes. For interesting ranges 
of Sm^* s and 5m^* e , the problematic neutrino evolution near Y e = 1/3 occurs well beyond 
the weak freeze-out radius and does not affect our treatment of the v* v^ v s and u* ^ u e 
transformations discussed previously. 

We have also ignored the possibility of a v e v^ v s resonance in the treatment of active- 
sterile neutrino transformation. This is because were this resonance to occur below the 
weak freeze-out radius, neutrino evolution through this resonance would be grossly non- 
adiabatic for the neutrino mass-squared difference (5m 2 es < 10~ 5 eV 2 ) adopted to solve the 



solar neutrino problem in our neutrino mixing scheme [see Eq. (12)]. 

To summarize the discussion in this subsection, we depict the squares of the effective 
neutrino masses as functions of density in Fig. 2. Here we assume that Y e > 1/3 and the 
order of the resonances is as shown. As Y e decreases, the effective mass track for v* will 
steepen, while that for v e will flatten out. For Y e < 1/3 this trend will be so extreme that 
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the order of the resonances will reverse. As noted above and discussed in Sec. IV, for our 
adopted neutrino mixing parameters this reversal does not occur below the weak freeze-out 
radius in a typical neutrino-heated outflow. 

C. Removal of the u e Flux and Evasion of the Alpha Effect 

Based on the discussion in Sec. IIIB, we can regard the the v* ^ v s and v* # v e reso- 
nances as filters of neutrino flavors/types when neutrino evolution through these resonances 
is adiabatic. Consider the evolution of an initial neutrino state 

\v{t = t )) = a e \u e ) +a s \u s ) + <v|^*) +a T *\v*). (24) 

At t > tfj,* s , i.e., after adiabatic propagation through the v* # v s resonance, the neutrino 
state becomes 

\v(t > t,,,*,)) « a' e \u e ) + a>;) + a^\u a ) + a' T *\v*), (25) 



where the primed coefficients differ from the corresponding unprimed ones in Eq. (|^ 
only by a phase. Symbolically, the function of the v* v^ v s resonance can be described as 
\u e ) — > \u e ), |z/ s ) — > 1 1/*), |z/*) — > \v s ), and \u*) — > \u*). At t > t M » e > £ M * S , i.e., after adiabatic 
propagation through the v* v^ v e resonance, the neutrino state becomes 

\u(t > V e )) « <>;) + o> e ) + <> s ) + aUK), (26) 

where the double-primed coefficients differ from the corresponding unprimed ones in Eq. 
(|4j) again only by a phase. Symbolically, the function of the v* v^ v e resonance can be 
described as \v e ) — > |z/*), \v s ) -> \v s ), \u*) — > \v e ), and |z/:) — > \u*). 

The probabilities for the initial neutrino state in Eq. ([24]) to evolve into the |z/ e ), |z/ s ), 
|f*), and \u*) states after adiabatic propagation through the z/* ^ v s and v* v^ f e resonances 
can be obtained from Eq. (|26|) as: 

^ e (t> Ve)~|(^k(t = to))f, (27) 
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P„.(t>Ve)« (KHt = t )) 



(28) 



^(* > Ve) « (^IO(^ek(i = to)) + (*/>;>(*/>(£ = to)) exp(z$) 



(v e \v(t = to)) + (z/>(t = t )) 



(29) 



and 



Pu r {t>t ll . 



(VrK)(VeW(t = t )) + (V r \v*)(u*\v(t = * )) 6Xp(l$) 



(z/ e |i/(t = t )) + (z/ T >(t = t )) 



(30) 



In Eqs. fl29|) and (|30|) the second approximation is obtained by averaging over the phase $. 
From Eqs. fl2"7l)- ([30~D we see that after adiabatic evolution through the u* v^ v s and 
v* v=^ z/ e resonances, (1) a z/ e emitted from the neutron star [i.e., \v{t = to)) = |f e )] would 
become 50% i/^ and 50% ^ T ; (2) a z/ M emitted from the neutron star [i.e., \u(t = to)) = 1^)] 
would become 50% u s , 25% u^, and 25% v T \ and (3) a i> r emitted from the neutron star 
[i.e., \v(t = to)) = \v T )] would become 50% u s , 25% v^, and 25% v T . Clearly, as no sterile 
neutrinos are emitted from the neutron star, the v e flux would be removed at those energies 
for which adiabatic evolution through the v* ^ u s and v* ^ v e resonances is completed. 
Therefore, if such neutrino evolution can be engineered to occur below the weak freeze-out 
radius and convert most of the v e emitted from the neutron star to other species before the 
electron fraction falls near or below 1/3, then we will have succeeded in disabling the alpha 
effect. 

IV. THE EVOLUTION OF THE ELECTRON FRACTION WITH 
ACTIVE-STERILE PLUS ACTIVE- ACTIVE NEUTRINO TRANSFORMATION 

Consider the Y e evolution of a mass/fluid element moving outward from the neutron star 
surface (i.e., the neutrino sphere). In the absence of neutrino flavor/type transformation, 
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this fluid element would be irradiated by the is e , z/ e , z/ M , i> M , i/ r , and v T fluxes emitted 
from the neutrino sphere throughout its progress toward ejection. (There should not be an 
appreciable flux of sterile neutrinos coming from the neutron star interior. Though a small 
sterile neutrino flux will not affect our conclusions, it is nonetheless interesting to note that 
for our assumed scheme of neutrino masses and mixings, matter effects inside the proto- 
neutron star should suppress strongly the production of sterile neutrinos.) Of course, the 
occurrence of the v* v^ v e resonance will cause the fluid element to experience deviation of 
its Y e evolution from the standard case with no neutrino flavor/type transformation. 

As discussed in Sec. IIIB, the v* v^ v s resonance occurs before (at higher temperature 
than) the v* v^ v e resonance for Y e > 1/3. Further, neutrinos with lower energies generally 
go through the two resonances before those with higher energies. Consequently, the two res- 
onances will sweep through the relevant neutrino energy distribution functions when viewed 
in the frame of a fluid element as it moves from the neutrino sphere at high temperature 
towards regions of lower temperature. The energy-position of the v* v^ v e resonance within 
the distribution functions, i.e., the resonance energy E^} ES , is related to temperature T 9 and 
radius r 6 through the resonance condition as 

0.4 



/Sin 
K RES « 29.4 MeV - M * e 



2.58 MeV 




10\ 3 

^100 



J 



_{8m?,* 



\ire 



\ 6eV 2 



'-'l00 r 6- (31) 



Y e + Y Ve -{Y^ + Y VT )/2\\ M NS 

When a fluid element reaches a given radius, we can regard all v e emitted from the 
neutrino sphere with energies less than the corresponding E^ ES at this radius as having 
changed to either v^ or v T . Note that no conversion into v e occurs as the \is*) component 
of the i/^ or v T emitted from the neutrino sphere with the relevant energies will have been 
converted into sterile neutrinos before reaching the v* v^ v e resonance. (Adiabatic neutrino 
evolution through the v* v^ v s and v* v^ v e resonances is assumed here.) The v e and 
V e capture rates in the fluid element, \ Uen and \ PeP , respectively, will both decrease with 
increasing radius due to the geometric dilution of the neutrino fluxes. However, as E^ ES 
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increases, conversion into v^ or v T will further decrease the v e flux in addition to the geometric 
dilution. Consequently, the rate A„ e „ for raising Y e will decrease more with increasing radius 
than the rate X PeP for lowering Y e . One would then expect Y e to be lowered progressively 
as the fluid element moves towards larger radius. If Y e < 1/3 is reached, the order of the 
v* t^ v s and v* v^ v e resonances will reverse (see Sec. IIIB) and further transformation of 
v e will cease. In principle this could be worrisome as the evasion of the alpha effect requires 
that the bulk of the v e flux over a substantial range of energies be transformed away prior 
to the point of alpha-particle formation. 

However, in practice the defeat of the alpha effect through removal of the v e flux will be 
achieved for the expected conditions in the neutrino-heated ejecta and plausible neutrino 
mixing parameters. This is because v e conversion coupled with the expected rapid expansion 
of the ejecta will cause the actual value of Y e at a given radius to differ significantly from 
the equilibrium value corresponding to the local v e and v e fluxes. This equilibrium value at 
radius r is defined as 

^ )5 u^r (32) 

As discussed above, in the presence of neutrino flavor/type transformation, the decrease of 
\ Uen due to conversion of v e into v^ or v T in addition to the geometric decrease common to 
both X UeTl and Ap eP tends to lower Y e E ^ at larger radius. Thus we have dY^/dr < 0. 

On the other hand, the Y e evolution of a fluid element prior to the point of alpha-particle 
formation is governed by 

dY (r) 
V{r) ^k = KAr) " [Kn{r) + A ^( r )]W, (33) 

where v (r) is the outflow velocity at radius r. To first order, the solution to Eq. (^) is 

As the fluid element moves to larger radius, the conversion "front" at E^ ES moves towards 
the high-energy end of the v e energy distribution function (shown schematically in Figure 



3), causing Y" e EQ to decrease. However, according to Eq. fl3~4|), this acts to increase Y e above 
yvEQ Further, a larger increase is obtained for a larger outflow velocity v corresponding 
to a more rapid expansion of the fluid element. Consequently, flavor evolution of v e with 
energies important for determining Y e can occur with Y e > 1/3 in the rapidly expanding 
neutrino-heated ejecta for a range of outflow conditions and neutrino mixing parameters. 

As indicated in Eq. (f31~l) , the resonance energy E^ ES at a given radius depends on the 
effective neutrino numbers relative to baryons Y Ue , Y Vfi , and Y Ut . These quantities determine 
the effects of the neutrino backgrounds on neutrino flavor/type evolution. Using Eq. @ 
and assuming adiabatic neutrino evolution, we find that for a radially travelling neutrino, 

10MeV\ [™ t /rn Wrn T /lOMeVY 



Y, 



A, 



^e,50 



.0MeV\ r°° /10M. 
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(35) 



where 
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0.723 



1.4M m \ 3 Slnr, 



Mi 
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(36) 



and 



fu e (E, 



El 



(37) 



TlF 2 ( V „ e )exp(E u jT Ue - Vl/e ) + l 

is the normalized z/ e energy distribution function. In Eqs. (|35|) and (|36|) 1^,50 and Lp e ^ 50 
are the neutrino energy luminosities in units of 10 50 ergs s _1 and R u $ is the neutrino-sphere 
radius in units of 10 6 cm. In Eq. ( |3"7| ) F 2 (r) Ve ) = J °° x 2 /[exp(x—r] Ue ) + l]dx is the second-order 
relativistic Fermi integral, and the two parameters T Ue and r] Ue can be specified by fitting 
the first two energy moments of the numerical v e energy spectrum fom supernova neutrino 
transport calculations. Expressions for Y v and Y Vt can be obtained in similar manner. 
For r > 1.5-Rj,, Eq. (Jffl) can be rewritten as 



A,. 



0.08 



T Q \ /1.4M Q ^ ' 



d2 q5 
rL u,6°100- 



(38) 



.107 V M NS 

Taking, for example, L Ve ~ L Pe ~ 10 50 erg s _1 , (E Ve ) ^ 11 MeV, and (E Pe ) ps 16 MeV, which 
are plausibly characteristic of neutrino emission at very late times, we have \Y V \ ~ 0.01 at 
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Tg ~ 10 for R v $ = 1 and S±oo = 0.7. Therefore, compared with typical values of Y e ~ 0.4, 
Y Ue , Y v , and Y^ can be neglected. Further, the effects of the neutrino backgrounds on 
neutrino flavor/type evolution are especially small for low values of Sioo [see Eq. (0)]. This 
is because important neutrino flavor/type evolution occurs over a relatively narrow range 
of temperatures and for a lower entropy, these temperatures correspond to larger radii [see 
Eq. (HD] where neutrino fluxes are smaller. 

In order to disable the alpha effect, we would like to have the bulk of the v e flux at 
E Ve < 20 MeV (see Fig. 3) transformed away before the point of alpha-particle formation 
(Tg ~ 10). This requires that the v e conversion front be at E^ ES « 20 MeV when the 
neutrino-heated ejecta has reached the radius corresponding to Tg ps 10. With a typical 
entropy of S'ioo = 0.7 and values of Y e ~ 0.4 3> Y Ve , Y v , Y Ut in the ejecta, we see from Eq. 
(pl|) that this requirement is indeed fulfilled for 5m 2 , e m 6 eV 2 . 

The typical evolution of Y e in the neutrino-heated ejecta with v* # v s plus v* # v e 
transformation can then be summarized as follows. At T 9 ~ 20, the electron fraction has a 
typical value of Y e ~ 0.4 in equilibrium with the v e and v e fluxes. As the ejecta moves to 
regions of lower temperature, the v e conversion front at E^ ES moves towards the high-energy 
end of the v e energy distribution function and removes the v e flux at E Ue < E^ ES . This 
suppresses the destruction of neutrons and lowers Y e somewhat below ~ 0.4. However, u e 
conversion coupled with the expansion of the ejecta causes rapid change in the instantaneous 
equilibrium Y e value, Y ECi , and maintains Y e > Y ECi . Consequently, conversion of v e with 
energies important for determining Y e is completed at Y e > 1/3. When the ejecta reaches the 
point of alpha-particle formation at Tg ~ 10, the u e conversion front is at E^ ES p» 20 MeV 
and the bulk of the v e flux has been removed. This then defeats the alpha effect which would 
drive Y e close to 0.5 if a significant flux of v e existed to destroy neutrons at T 9 < 10. 

Therefore, our scenario of v* v^ v s plus v* # v e transformation ensures that 1/3 < 
Y e < 0.4 can be obtained and maintained before the onset of rapid neutron capture in 
the neutrino-heated ejecta. While this scenario does not lead to very low values of Y e , it 
nevertheless guarantees a range of Y e that will lead to a successful r-process when combined 
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with the appropriate entropy S and the dynamic expansion timescale tdyn in the ejecta. 

V. CONCLUSIONS 

The alpha effect is the single biggest impediment to obtaining the necessary conditions for 
r-process nucleosynthesis in neutrino-heated ejecta from supernovae. In fact, this effect may 
also be an obstacle for r-process nucleosynthesis in neutrino-heated ejecta from neutron-star 
mergers. The central problem is that the v e flux causing ejection of baryons from deep in 
the gravitational potential well of the neutron star will destroy neutrons via v e + n —>■ p + e~ 
in the regions of alpha-particle formation. This destruction of neutrons then renders the 
subsequent neutron capture process incapable of producing the heavy r-process elements. We 
have suggested a scheme of neutrino masses and mixings which could solve this conundrum by 
removing the v e flux through matter-enhanced v^ v^ v s plus v li>T v^ u e transformation above 
the regions of efficient neutrino heating but below the regions of alpha-particle formation. 
This scheme which rescues r-process nucleosynthesis in neutrino-heated ejecta originally was 
not constructed for this purpose. Rather, it was designed to explain simultaneously the solar 
neutrino data and the anomalous atmospheric v^jv,, ratio and to allow for a hot component 
of dark matter. Subsequently, it also explained the LSND signal. In fact, this scheme with a 
maximally-mixed v il -v r doublet split from a lower-mass v e -v a doublet by (<5m 2 )LSND ^ leV 2 
may be the only one which can escape elimination by recent interpretations of the SuperK 
atmospheric neutrino data |7,25|1 and by Big Bang Nucleosynthesis considerations ||34| . 



While we do not know if there is a way other than invoking neutrino mixing to circumvent 
the alpha effect in neutrino-heated ejecta, attempts at "astrophysical" dodge of this effect 
made so far (see, e.g., Ref. |llj) seem to be finely tuned at best. Absent a non-neutrino- 
physics escape from the alpha effect, we could draw interesting inferences about Galactic 
chemical evolution were future experiments to reveal a neutrino mass scheme other than 
the one that aids the r-process. In that case, for example, we could be forced to re-think 
the origin of the bulk of the r-process material in the Galaxy. In turn, this could have 
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consequences for our understanding of the rates and physics of, e.g., neutron-star mergers. 
Or we could be forced to re-think electromagnetic ejection of material from supernovae. 

In addition to the r-process connection, our scheme of neutrino masses and mixings has 
interesting consequences for detection of neutrinos from future Galactic supernovae. As the 
solar neutrino problem is solved by v e v^ v s mixing in this scheme, sterile neutrinos produced 
by the v^ T v^ v a transformation near the neutron star will be converted into v e . Therefore, 
at late times of the supernova process, there may be a significant increase in the average 
v e energy with no accompanying increase in the average v e energy. Furthermore, it has 
been shown that neutrino flavor transformation has important effects on the dynamics of 
supernova explosion p9|-p2|. Conceivably, neutrino mixing also affects the nucleosynthesis 



discussed in Ref . |35j , which occurs shortly after the supernova explosion. Currently, we are 
investigating the effects of our scheme of neutrino masses and mixings at these early times. 
From this work, one conclusion specific to particle physics is that evasion of the alpha 
effect and, hence, robust production of the r-process elements in neutrino-heated ejecta seem 
to require at least one light sterile neutrino species. We have shown here that one neutrino 
mass scheme to disable the alpha effect has a maximally-mixed v iX -v T doublet split from 



a lower-mass v e -v s doublet. However, another scheme |33|| , which has three light, nearly- 
degenerate active neutrinos and a heavier sterile neutrino species, has also been suggested to 
do the same. Though this latter scheme has the attractive feature of a significant increase 
in the neutron excess, considerations that do not concern the r-process argue against it. 
Finally, if we adopt the neutrino mass scheme of this paper, then the splitting between 
the doublets must be 5m? > 1 eV 2 in order to have the beneficial effects on the r-process. 
Note that this splitting is within the LSND range and most likely in the upper sector of it. 
Of course, these conclusions presuppose that some of the r-process material in the Galaxy 
originated in neutrino-heated ejecta (from either supernovae or neutron-star mergers) and 
that there is no conventional astrophysical fix for the alpha effect. 

Our conclusions will be tested by neutrino experiments as well as and astrophysical 
observations on r-process nucleosynthesis. In any case, it is both surprising and tantalizing 
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that mixings among active and sterile neutrino species with small masses can have such 
profound effects on the physics of astrophysical objects and the synthesis of the heaviest 
elements. 
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FIGURES 
FIG. 1. The neutrino mass scheme discussed in this paper. A doublet of 

(near)-maximally-mixed v^ and v T neutrinos with mass-squared difference <5m 2 ~ 10~ 2 eV 2 is 

split from a doublet of lower- mass v e and v s by a mass-squared difference 5m|} oublets . 

The lower-mass doublet can be arranged (lower right) with maximal (or near maximal) 

v e v^ v s vacuum mixing to give a "just so" solar neutrino solution, which would require 

5m? es ~ 10~ 10 eV 2 . Alternatively (lower left), the sterile neutrino could be heavier than the 

v e with 5m 2 es < 1CT 5 eV 2 to give a matter-enhanced (MSW) solution to the solar neutrino 

problem. We assume that <5m 2 oublets is much larger than the mass-squared splittings within 

the doublets. 

FIG. 2. Cartoon of the instantaneous neutrino mass levels (effective mass-squared m 2 ff ) as 
functions of matter density p for Y e > 1/3. 

FIG. 3. Example v e and iv T energy distribution functions. Here we take these functions to be 
of the form in Eq. ( p7| ) with T Ue = 2.75 MeV, r] Ve = 3 (corresponding to (E Ue ) = HMeV) and 
T v iT = 6.76 MeV, r) v = 3 (corresponding to (E u T ) = 27 MeV). The resonance energy E^- ES for 
neutrino flavor /type transformation sweeps from low to high energy through the neutrino energy 
distribution functions as a fluid element moves away from the neutron star. 



28 



TABLES 
TABLE I. Weak potentials derived from neutrino forward scattering on electron/nucleon [A) 

and neutrino (B) backgrounds for various channels of neutrino transformation. The corresponding 

antineutrino transformation channels have opposite signs for A and B. Here N A is Avogadro's 

number, G F is the Fermi constant, p is the matter density, Y e is the electron fraction, and Y Ue , 

Y v , and Y Ut are the effective neutrino numbers for the corresponding species relative to baryons. 



Channel 


A 




B 


v e # v s 


(3V2/2)G FP N A (Y e - 


-1/3) 


V2G FP N A (2Y Ue + Y^ + Y UT ) 


»H ^ V s 


(V2/2)G FP N A (Y e - 


-1) 


V2G F pN A (Y„ e + 2Y^ + Y UT ) 


V T ^ v s 


(V2/2)G FP N A (Y e - 


-1) 


y/2G F pN A (Y Ve +Y Vll + 2Y UT ) 


"H # V e 


s/2G F pN A Y e 




V2G FP N A (Y Ue - yj 


V T # v e 


V2G F pN A Y e 




V2G FP N A (Y Ue -Y UT ) 



Vfj, v- V T 



V2G FP N A (Y^-Y Ur ) 
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<Sm 2 _ r ~ 1(T 2 eV 2 



V 



V 



V T 



8m 



doublets 



v s 



- ve 



Sml e < 10" 5 eV 2 6m*_ s ~ 10" 10 eV 2 




Fast expansion case: 



followed by ^//,t ^ v e 



In the fast expansion case, y e lags neutrino transformation. 
Therefore, the order of the level crossings does not change 
until all of the electron neutrinos are transformed. 



Net Result: Probability v li "^ v li 
Probability ^ji ~^ v T 
Probability ^ji ~^ v s 



1/4 
1/4 

1/2 



Probability ^/x,r - *■ v e =>■ (2nd resonance) 



/; 



V 



0.08 



0.06 



0.04 



0.02 




Vt,vh 



E v {MeV) 



